Search for: All records

ABSTRACT Fault zones exhibit geometrical complexity and are often surrounded by multiscale fracture networks within their damage zones, potentially influencing rupture dynamics and near-field ground motions. In this study, we investigate the ground-motion characteristics of cascading ruptures across damage zone fracture networks of moderate-size earthquakes (Mw 5.5–6.0) using high-resolution 3D dynamic rupture simulations. Our models feature a listric normal fault surrounded by more than 800 fractures, emulating a major fault and its associated damage zone. We analyze three cases: a cascading rupture propagating within the fracture network (Mw 5.5), a non-cascading main-fault rupture with off-fault fracture slip (Mw 6.0), and a main-fault rupture without a fracture network (Mw 6.0). Cascading ruptures within the fracture network produce distinct ground-motion signatures with enriched high-frequency content, arising from simultaneous slip of multiple fractures and parts of the main fault, resembling source coda-wave-like signatures. This case shows elevated near-field characteristic frequency (fc) and stress drop, approximately an order of magnitude higher than the estimation directly on the fault of the dynamic rupture simulation. The inferred fc of the modeled vertical ground-motion components reflects the complexity of the radiation pattern and rupture directivity of fracture-network cascading earthquakes. We show that this is consistent with observations of strong azimuthal dependence of corner frequency in the 2009–2016 central Apennines, Italy, earthquake, sequence. Simulated ground motions from fracture-network cascading ruptures also show pronounced azimuthal variations in peak ground acceleration (PGA), peak ground velocity, and pseudospectral acceleration, with average PGA nearly double that of the non-cascading cases. Cascading ruptures radiate high-frequency seismic energy, yield nontypical ground-motion characteristics including coda-wave-like signatures, and may result in a significantly higher seismologically inferred stress drop and PGA. Such outcomes emphasize the critical role of fault-zone complexity in affecting rupture dynamics and seismic radiation and have important implications for physics-based seismic hazard assessment. 

Abstract Seismic faults are surrounded by damaged rocks with reduced rigidity and enhanced attenuation. These damaged fault zone structures can amplify seismic waves and affect earthquake dynamics, yet they are typically omitted in physics‐based regional ground motion simulations. We report on the significant effects of a shallow, flower‐shaped fault zone in foreshock‐mainshock 3D dynamic rupture models of the 2019 Ridgecrest earthquake sequence. We find that the fault zone structure both amplifies and reduces ground motions not only locally but at distances exceeding 100 km. This impact on ground motions is frequency‐ and magnitude‐dependent, particularly affecting higher frequency ground motions from the foreshock because its corner frequency is closer to the fault zone's fundamental eigenfrequency. Within the fault zone, the shallow transition to a velocity‐strengthening frictional regime leads to a depth‐dependent peak slip rate increase of up to 70% and confines fault zone‐induced supershear transitions mostly to the fault zone's velocity‐weakening roots. However, the interplay of fault zone waves, free surface reflections, and rupture directivity can generate localized supershear rupture, even in narrow velocity‐strengthening regions, which are typically thought to inhibit supershear rupture. This study demonstrates that shallow fault zone structures may significantly affect intermediate‐ and far‐field ground motions and cause localized supershear rupture penetrating into velocity‐strengthening regions, with important implications for seismic hazard assessment. 

Several regularly recurring moderate-size earthquakes motivated dense instrumentation of the Parkfield section of the San Andreas fault, providing an invaluable near-fault observatory. We present a seismo-geodetic dynamic inversion of the 2004 Parkfield earthquake, which illuminates the interlinked complexity of faulting across time scales. Using fast-velocity-weakening rate-and-state friction, we jointly model 3D coseismic dynamic rupture and the 90-day evolution of postseismic slip. We utilize a parallel tempering Markov chain Monte Carlo approach to solve this non-linear high-dimensional inverse problem, constraining spatially varying prestress and fault friction parameters by 30 strong motion and 12 GPS stations. From visiting >2 million models, we discern complex coseismic rupture dynamics that transition from a strongly radiating pulse-like phase to a mildly radiating crack-like phase. Both coseismic phases are separated by a shallow strength barrier that nearly arrests rupture and leads to a gap in the afterslip. Coseismic rupture termination involves distinct arrest mechanisms that imprint on afterslip kinematics. A backward propagating afterslip front may drive delayed aftershock activity above the hypocenter. Analysis of the 10,500 best-fitting models uncovers local correlations between prestress levels and the reference friction coefficient, alongside an anticorrelation between prestress and rate-state parameters b−a. We find that a complex, fault-local interplay of dynamic parameters determines the nucleation, propagation, and arrest of both, co- and postseismic faulting. This study demonstrates the potential of inverse physics-based modeling to reveal novel insights and detailed characterizations of well-recorded earthquakes. 

Abstract With the rise of data volume and computing power, seismological research requires more advanced skills in data processing, numerical methods, and parallel computing. We present the experience of conducting training workshops in various forms of delivery to support the adoption of large-scale high-performance computing (HPC) and cloud computing, advancing seismological research. The seismological foci were on earthquake source parameter estimation in catalogs, forward and adjoint wavefield simulations in 2D and 3D at local, regional, and global scales, earthquake dynamics, ambient noise seismology, and machine learning. This contribution describes the series of workshops delivered as part of research projects, the learning outcomes for participants, and lessons learned by the instructors. Our curriculum was grounded on open and reproducible science, large-scale scientific computing and data mining, and computing infrastructure (access and usage) for HPC and the cloud. We also describe the types of teaching materials that have proven beneficial to the instruction and the sustainability of the program. We propose guidelines to deliver future workshops on these topics. 

Abstract Dynamic rupture simulations generate synthetic waveforms that account for nonlinear source and path complexity. Here, we analyze millions of spatially dense waveforms from 3D dynamic rupture simulations in a novel way to illuminate the spectral fingerprints of earthquake physics. We define a Brune-type equivalent near-field corner frequency (fc) to analyze the spatial variability of ground-motion spectra and unravel their link to source complexity. We first investigate a simple 3D strike-slip setup, including an asperity and a barrier, and illustrate basic relations between source properties and fc variations. Next, we analyze >13,000,000 synthetic near-field strong-motion waveforms generated in three high-resolution dynamic rupture simulations of real earthquakes, the 2019 Mw 7.1 Ridgecrest mainshock, the Mw 6.4 Searles Valley foreshock, and the 1992 Mw 7.3 Landers earthquake. All scenarios consider 3D fault geometries, topography, off-fault plasticity, viscoelastic attenuation, and 3D velocity structure and resolve frequencies up to 1–2 Hz. Our analysis reveals pronounced and localized patterns of elevated fc, specifically in the vertical components. We validate such fc variability with observed near-fault spectra. Using isochrone analysis, we identify the complex dynamic mechanisms that explain rays of elevated fc and cause unexpectedly impulsive, localized, vertical ground motions. Although the high vertical frequencies are also associated with path effects, rupture directivity, and coalescence of multiple rupture fronts, we show that they are dominantly caused by rake-rotated surface-breaking rupture fronts that decelerate due to fault heterogeneities or geometric complexity. Our findings highlight the potential of spatially dense ground-motion observations to further our understanding of earthquake physics directly from near-field data. Observed near-field fc variability may inform on directivity, surface rupture, and slip segmentation. Physics-based models can identify “what to look for,” for example, in the potentially vast amount of near-field large array or distributed acoustic sensing data. 

Abstract Several regularly recurring moderate‐size earthquakes motivated dense instrumentation of the Parkfield section of the San Andreas fault (SAF), providing an invaluable near‐fault observatory. We present a seismo‐geodetic dynamic inversion of the 2004 Parkfield earthquake, which illuminates the interlinked complexity of faulting across time scales. Using fast‐velocity‐weakening rate‐and‐state friction, we jointly model coseismic dynamic rupture and the 90‐day evolution of postseismic slip in a 3D domain. We utilize a parallel tempering Markov chain Monte Carlo approach to solve this non‐linear high‐dimensional inverse problem, constraining spatially varying prestress and fault friction parameters by 30 strong motion and 12 GPS stations. From visiting 2 million models, we discern complex coseismic rupture dynamics that transition from a strongly radiating pulse‐like phase to a mildly radiating crack‐like phase. Both coseismic phases are separated by a shallow strength barrier that nearly arrests rupture and leads to a gap in the afterslip, reflecting the geologic heterogeneity along this segment of the SAF. Coseismic rupture termination involves distinct arrest mechanisms that imprint on afterslip kinematics. A backward propagating afterslip front may drive delayed aftershock activity above the hypocenter. Trade‐off analysis of the 10,500 best‐fitting models uncovers local correlations between prestress levels and the reference friction coefficient, alongside an anticorrelation between prestress and rate‐state parameters . We find that a complex, fault‐local interplay of dynamic parameters determines the nucleation, propagation, and arrest of both, co‐ and postseismic faulting. This study demonstrates the potential of inverse physics‐based modeling to reveal novel insights and detailed characterizations of well‐recorded earthquakes.